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Counting spanning trees in self-similar networks by evaluating determinants 
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Spanning trees are relevant to various aspects of networks. Generally, the number of spanning 
trees in a network can be obtained by computing a related determinant of the Laplacian matrix of the 
network. However, for a large generic network, evaluating the relevant determinant is computation- 
ally intractable. In this paper, we develop a fairly generic technique for computing determinants 
corresponding to self-similar networks, thereby providing a method to determine the numbers of 
spanning trees in networks exhibiting self-similarity. We describe the computation process with a 
family of networks, called [x, y)-flowers, which display rich behavior as observed in a large variety 
of real systems. The enumeration of spanning trees is based on the relationship between the de- 
terminants of submatrices of the Laplacian matrix corresponding to the (x, t/)-flowers at different 
generations and is devoid of the direct laborious computation of determinants. Using the proposed 
method, we derive analytically the exact number of spanning trees in the (x, y)-flowers, on the basis 
of which we also obtain the entropies of the spanning trees in these networks. Moreover, to illus- 
trate the universality of our technique, we apply it to some other self-similar networks with distinct 
degree distributions, and obtain explicit solutions to the numbers of spanning trees and their en- 
tropies. Finally, we compare our results for networks with the same average degree but different 
structural properties, such as degree distribution and fractal dimension, and uncover the effect of 
these topological features on the number of spanning trees. 

PACS numbers: 89.75.Hc, 05.50.+q, 05.20.-y, 04.20.Jb 
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I. INTRODUCTION 

Counting spanning trees in networks is a fundamental 
issue in both theory and applications, which has recently 
attracted increasing attention from mathematics [IHj], 
physics [EI-fToj , and other fields [ll[ • As a network invari- 
ant, span ning tree is a crucial measure of the network 
reliability |L^ - tl4| . The notion is also closely related to 
various aspects of networks. For example, the number 
of spanning trees of a connected network is exactly the 
number of recurrent configurations of the Abelian sand- 
pile model on the network pl| [l6| , which is equivalent to 
the chip-firing game [lTj under certain restrictions and is 
a paradigm for self-organized criticality [l8[. As another 
example, spanning trees offer useful insights into under- 
standing the origin and mechanism of fractality in fractal 
scale- free networks [19j . In addition, spanning trees are 
relevant to other interesting problems on networks, such 
as transpor t l20ll. e lectrical networks [2l| . unbiased ran- 
dom walks [22L |23|. loop-erased random walks [Hj], and 
so on. 

In view of their relevance to diverse aspects of net- 
works and a wide range of applications [25| . spanning 
trees in networks have become a focus of some recent 
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research [26H32J. Particularly, in the physics literature 
a lot of efforts have been devoted to enumerating span- 
ning trees in specific networks by using different tech- 
niques according to their special structures. Examples 
include regular lattices [1-0, S, Sierpinski gaskets d, , 
Erdsd-Renyi random graphs [4J, the pseudofractal scale- 
free web [33j], fractal scale- free networks [Hj], and so 
on. These works have provided some effective methods 
for determining spanning trees in special networks, and 
have revealed some nontrivial effects of network struc- 
tural properties on spanning trees. 

Most existing methods for counting the number of 
spanning trees are applicable only to particular networks. 
Due to the complexity and diversity of networks, a gen- 
eral approach for counting spanning trees of a generic 
network is not available, with the exception of the Kirch- 
hoff Matrix Tree Theorem [35| . This well-known theorem 
provides a universal algorithm for determining the num- 
ber of spanning trees of an arbitrary connected graph, 
in terms of a determinant [36]. However, the compu- 
tational complexity for evaluating the determinant of a 
general network is very high |37| . For large networks, it 
is difficult and even impossible to obtain their numbers of 
spanning trees. For this reason, seeking an efficient and 
somewhat general technique for calculating the determi- 
nant of a network is a matter of exceptional importance. 

On the other hand, it has been observed that by na- 
ture most real networks are fractal [38| as well as scale- 
free [39l ] . Both striking properties have significant in- 
fluences on other structural properties, e.g., average dis- 
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tance (40l - |42l | and degree correlations [H}. In particular, 
existing results have shown that in some special networks, 
the scale- free property may decrease the number of span- 
ning trees [33[ , while fractality can increase the number of 
spanning trees [34j. Thus, some interesting questions rise 
naturally: Can the power-law behavior alone determine 
the number of spanning trees in a scale-free network? 
What is the relation between the fractal dimension and 
the spanning trees in a fractal scale-free network? Both 
questions remain unsolved today. 

In this paper, we present an approach for explicitly 
determining the number of spanning trees by evaluat- 
ing a related determinant in a network, which is ex- 
pected to hold true for a broad class of self-similar graphs. 
Our method is based on the recursion relations of de- 
terminants of sub-matrices of the Laplacian matrix as- 
sociated with the network at different iterations, which 
avoids the laborious computation of the largest deter- 
minant as needed by classic algebraic schemes. To the 
best of our knowledge, this constitutes the first work in 
exploring fast evaluation of the determinant of a com- 
plex network, especially for networks with scale- free [3!| , 
small- world [44j> and fractal [38| properties. 

To demonstrate the computation process of our tech- 
nique, we apply it to find the number of spanning 
trees in the family of self-similar networks, called (x, y)- 
fiowers [ill HH , which display some remarkable proper- 
ties observed in real-life networks [47l - l49| and include the 
extensively studied pseudofractal scale-free web [H(3, [H| 
and fractal hierarchical lattices [52l[55j as their limiting 
cases. Using the proposed technique, we derive the exact 
number of spanning trees in the whole family of [x, y)- 
fiowers, based on which we obtain their entropies of span- 
ning trees. In addition, in order to exhibit the generality 
of our method, we apply it to some other self-similar net- 
works and obtain closed-form solutions for the numbers 
of spanning trees in such networks. Finally, we present a 
detailed analysis of the obtained results, and show that 
power-law distribution alone do not suffice to character- 
ize the numbers of spanning trees in scale-free networks, 
and that the increasing fractal dimensions of fractal scale- 
free networks may lead to decrements of the numbers of 
spanning trees. 



II. MODEL AND PROPERTIES OF 

(x,y)-FLOWERS 

In this section, we introduce the construction and 
structural properties of the (x, y) -flowers 0, [46jj, which 
are built in an iterative way. Let F n (x,y) (n > 0) de- 
note the (x, y)-flowers after n generations of evolution. 
In what follows it is assumed that x < y and y > 1 
without loss of generality. The (x, y)-flowers are con- 
structed as follows. For n — 0, F (x,y) is an edge con- 
necting two nodes, called the initial nodes hereafter. For 
n > 1, F n (x,y) is derived from F n ^\(x,y) through re- 
placing each exiting edge in F n -i(x,y) by two parallel 



paths of lengths x and y, see Fig. [T] For illustration, in 
Figs. [2] and G2 we present the growing processes of two 
particular networks: (l,3)-flower and (2, 2)-flower. Ac- 
cording to the generating algorithm for the (x, y)-flowers, 
it is easy to see that the number of edges in F n (x, y) is 

M n = {x + y) n . (1) 




FIG. 1: (Color online) Iterative construction method of the 
(x, y)-flowers. Each edge is replaced by two parallel paths 
with lengths x (x > 1) and y (y > x and y > 1) on the right- 
hand side of the arrow. For x = 1, a pair of old nodes directly 
connected by an old edge generates y — 1 new nodes. All of 
these y — 1 new nodes and the two old nodes form a red path 
of length y; while the old edge is remained. For x > 1, each 
old edge connecting two old nodes is removed and replaced 
by two paths with the two old nodes as the ends of the two 
paths: the y — 1 red nodes and the two old nodes form a path 
of length y, while the x — 1 green nodes and the two old nodes 
constitute another path of length x. 

Alternatively, the (x, y)-flowers can be constructed 
by using the following equivalent algorithm [45j . which 
highlights the self-similarity of the networks. From 
F n -i(x, y), the next generation F n (x,y) is obtained by 
joining x + y copies of F n -i(x.y) at the two initial 
nodes of each replica. Figure [4] provides an illustra- 
tion for the second construction approach for the (1,3)- 
flower. According to this construction, the number of 
nodes N n in F n (x, y) satisfies the recursive relation N n — 




n=3 



FIG. 2: (Color online) Illustration of the growing process for 
the (1, 3)-flower. 
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n=0 




family of F n (x,y), all networks with the same x + y have 
an identical degree sequence, hence an identical degree 
distribution. 

The (x, ;/ )-no wers exhibit rich behavior in their archi- 
tecture [45, 46] . They obey a power-law degree distribu- 
tion P(k) ~ fc~ 7 with the exponent 7 = l + ln(x + j/)/ln2 
lying in the interval [l + ln3/ ln2, 00). For x = 1, the net- 
works are small- world but non- fractal; while for x > 1, 
they are "large- world" and fractal with the fractal dimen- 
sion df = ln(x + y)/ hix. These topological properties are 
not shared by any other modeling networks, which make 
the (x, j/)-flowers unique within the family of networks. 




n=3 

FIG. 3: (Color online) Sketch of the iterative process for the 
(2, 2)-flower. 



(x+y)N n —i — (x+y), which together with the initial con- 
dition iVi = x + y yields 



N n = X -±^(x + yT + - X + y 
1 x 



x + y — 1 " x + y — 1 

Thus, the average degree of the (a;, y)-flowers is 

2M n 2(x + y~l)(x + yr- 1 
< fc >™ = Tr- = 



(2) 



N n {x + y - 2){x + y) n - 1 + 1 ' 
which approaches 2(x + y — l)/(x + y — 2) as n — > 00. 



(3) 




I Fg(l;3 J I F*(l;3 ] ( |F: 




F n+ 1 (1 ;3 ) 



FIG. 4: (Color online) Illustration of the alternative construc- 
tion algorithm for the (1, 3)-flower. F n +i(l,3) is obtained 
through amalgamating four copies of F n (l,3), denoted by 
F%(1,3) (rj = 1,2,3,4), by merging four pairs of the initial 
nodes belonging to different F%(1,3). Black nodes represent 
the initial nodes of F%(1, 3) or of F n +i(l, 3). 

The deterministic construction also allows to deter- 
mine the full degree distribution of the (x, y)-flowers. In 
F n (x, y), the degree spectra of all nodes are discrete, i.e., 
nodes have only degree k in the form of k — 2 m (m = 
1,2, ... , n). Let N n (m) be the number of nodes with de- 
gree 2 m in F n (x,y). Then, N n (m) = {x+y-2)(x+y) n ~ m 
for m < n, and N n (m) = x + y for m = n. Thus, in the 



III. SPANNING TREES IN (i,y)-FLOWERS 

After introducing the construction and some character- 
istics of the (x, y)-flowers, we now attempt to determine 
analytically the number of spanning trees in F n (x,y), 
denoted by NstM. According to the Matrix Tree Theo- 
rem [35], [3y|, NsT{n) can be evaluated by computing the 
determinant of a submatrix of the Laplacian matrix cor- 
responding to F n (x,y), which is obtained by removing 
a row and a column corresponding to an arbitrary node 
of F n (x,y). The elements of the Laplacian matrix asso- 
ciated with F n (x,y), represented by F n , are defined as 
follows: the non-diagonal entry (i ^ j) is —1 (or 0) if 
nodes i and j are (or not) directly connected by an edge, 
while the diagonal entry fa equals the degree of node i. 

To find Nst(ti), we denote by G n the submatrix of 
F n obtained by removing from F n the row and column 
corresponding to an initial node of F n (x,y). Then, we 
have 



AT ST (n) = det(G„). 



(4) 



Our aim is to determine det(G„). Although for a general 
network of large size, the evaluation of its determinant 
is very hard and even impossible, below we show that 
for F n (x, y) and other self-similar networks, we can eval- 
uate the associated determinants and derive closed-form 
expressions for the number of their spanning trees. 

First, we define another submatrix of F n , which is ob- 
tained by deleting the two rows and columns from F n 
corresponding to the two initial nodes, and denote it by 
H n . In what follows, we show the recursive relations 
for det(G n ) and det(iT„), thereby, expressing det(G„+i) 
and det(i? n _|_i) in terms of det(G n ) and det(-ff„), based 
on which we offer explicit solutions for these two quan- 
tities. For simplicity, let /„, g n , and h n represent 
det(F„), det(G„), and det(H n ), respectively. Obviously, 
/„ = det(.F n ) = for all n, since the sum of each row in 
F n is zero. 

It should be noted that the technique we develop 
for evaluating the determinants can be adapted to a 
wide class of self-similar networks other than the (x, y)- 
flowers. Here, we illustrate the computation process with 
the (a;, y)-flowers. Moreover, we only consider in detail 
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two particular and representative networks in the (x, y)- 
flower family, namely, (1,3) -flower and (2, 2) -flower. For 
the general (x, y)-flowers we can treat them analogously 
although the computation is somewhat complex. In so 
doing, for notational simplicity, we only provide the re- 
cursion relations governing g n and h n for the general 
(x, y)-flowers, but omit their detail computations. 



A. The (l,3)-flower 

According to the second construction algorithm of the 
networks, see Fig. 01 the Laplacian matrix of F„ + i(l,3), 
denoted by F n +i in the case without confusion, can be 
expressed in terms of H n as 
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In Eq. ([5]), the first two rows and two columns correspond 
to the two initial nodes of 3); d n is the degree of 

either of initial nodes at generation n; p n and q n are two 
row vectors with order N n — 2 = §(4" — 1) describing re- 
spectively the interaction between the two initial nodes 
and the other nodes in different replicas of F n (l, 3) that 
are joined together forming F n +i{l,3); and the super- 
script T of a vector represents its transpose. The third 
and forth rows and columns describe the connection re- 
lation of the two nodes created at n = 1 (each of which 
has the same degree 2d n as the two initial nodes) and 
other nodes in F n+ i(l, 3). 

Based on Eq. ((SJ), we can also represent G„+i and H n+ i 
in terms of H n , as 
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2d 



H 



n+l 



-1 
1 2d n 



Qn Pn 
q n 



T 


T 


qZ 


T 


pZ 


Qn 




pZ 



H n O 

O H r 
O 
O 



\ 

Pn 

o 
o 

O H n O 
O O H n ) 



o 
o 



(7) 



In the case of n = 2, for example, Eqs. 
become respectively 
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With the expressions for G n +i and H n +i given in 
Eqs. (O and ([7]), we can evaluate the determinants of 
G n and H n , namely, g n and h n , on the basis of their 
recursive relations that will be derived. We first show a 
detailed computation process for g n , while the evaluation 
for h n can be dealt with in an analogous way. 

Note that the first row of G n +i can be viewed as the 
sum, i.e., a linear combination with scalars being 1, of 
two vectors; that is, 

( 2d n -1 9 „J)„0 0) 
= ( d n q n ) + ( d n -1 p n ) . 

Analogously, the second and third rows can also be re- 
garded separately as the sum of two vectors, namely, 

( -1 2d„ -1 q n Pn 0) 
= ( -1 d n q n ) + ( d n -1 Q p n ) , 

and 

(0 -1 2d„ q n Pn ) 
= ( -1 d n q n ) + ( d n pn ) . 

Thus, according to the properties of determinants, g n +i 
can be decomposed as the sum of eight determinants, as 
follows: 
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Therefore, the evaluation of g n +i is reduced to finding 

the values of the eight determinants on the right-hand 

side of Eq. @, which are denoted sequentially as 

(1 < % < 8). Using some elementary matrix operations, 

(1) - o (2) - n {4) - o (8) and 
— 9n+l — 9n + l — 9n+l ana 

(7) 

g n _l_ v Having g n +i m terms 
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of the quantities of g^+ii the next step is to explicitly 
determine these quantities. In view of the equivalence of 
some g n l) +l , we only compute g n 2 \ x and g^. 

(2) (3) 

For g n i-\ and utilizing some elementary matrix 

operations, we can easily obtain 
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Plugging these expressions into Eq. (|8]), we have 
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derive the recursion 
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Next, we derive the recursion relation for 
Through a computational process similar to that for 
(?„+i, we obtain the following relation: 
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Having obtained the above recursive expressions for 9n+i an d h n +i, we are now in a position to explicitly de- 
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tcrminc g n . For this purpose, we examine a new quantity 
9n/hn, which obviously satisfies the following recursive 
relation: 

h n+ i 3 h n 

Considering the initial condition g\jh\ = 4/3, the above 
equation is solved by induction to yield 



respectively. In Eqs. (fTTj) and (IT21) . r„ and s n are two vec- 
tors describing the connection relation between an initial 
node and the other nodes. 

Based on the above two relations, we obtain the follow- 
ing relations for the determinants corresponding to G n +i 
and H n+ i by using some elementary matrix operations: 



9n 

h n 



Hence, 



(9) 



Substituting Eq. ([§]) into Eq. flSJ) and using gi — 4, we 
can solve Eq. ([8]) to obtain the exact number of spanning 
trees: 

N ST (n) = g n = 3 C4"-3n-l)/9 4 (2x4»+3„-2)/9 > (1Q) 

It is easy to express Nst(ti) in terms of the network 
size N n , so as to obtain the relation governing the two 
quantities. Recalling that for (1, 3)— flower, N n = (2 X 
4" + 4)/3, hence we have 4" = (3N n - 4)/2 and n = 
[ln(3iV„ - 4)/ In 2 - l]/2. These relations allow us to 
recast NsT{n) as a function of N n : 

N ST (n) = 3[ JV "- ln (3A r n-4)/ln2-l]/6 4 [2AT„+ln(3JV„-4)/ln2-5]/6 

which indicates that Nst(ji) approximately increases ex- 
ponentially in N n . 

Then we can define the entropy of spanning trees for 
the (1, 3)— flower as the limiting value 
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= -(41n2+ln3) « 0.6452. 
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B. The (2,2)-flower 

We proceed to determine the number of spanning trees 
in the (2, 2)-flower. In the computation process, we use 
the same notations as those defined for the (1, 3)-flower. 
Similarly to the (1, 3)-flower, we can derive the recursion 
relations for G„+i and H n+ i, as 
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and 
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9n+i = ^K{g n f + 4(/i„) 2 /„g„ = Ah n (g n f, 



and 



±(h n f{g n f 



These two recursion relations, coupled with the initial 
conditions g\ = 4 and hi = 4, yield 



N ST (n)=g n = 2i^- 1 \ 



(13) 



which can be further expressed in terms of the network 
size N n as 



N ST (n) = 2 N -- 2 . 

Thus, for the (2, 2)-flower, the number of spanning trees 
also scales exponentially with the network size. And its 
entropy of spanning trees is 



E, 



(2,2) -flower 



= lim 



In Nst (n) 

Nn 



= In 2 w 0.6931 . 



a value larger than that corresponding to the (1,3)- 
flower. 



C. The (x, y)-flowers 

The above process for computing spanning trees in the 
(l,3)-flower and the (2,2)-flower can be used to deter- 
mine the number of spanning trees in the whole family 
of (x, y)-flowers. Since for the general (x, y)-flowers, the 
recursive relations for G„+i and H n+ \ arc lengthy and 
the computation process for g n +i and h n +\ is more com- 
plex, we only provide the recursive relations for g n +i and 
h n +i, leaving out the derivation detail, which is similar to 
those for the two particular cases considered above. The 
recursive relations of g n +i and h n+1 for the (x, y)-flowers 
are given by 



x+y—l 



and 



9n+i = (x + y)h n {g n ) 



h n+1 =xy(h n ) 2 (g n ) x+ y- 2 . 



(14) 



(15) 



Considering the initial conditions g\ — x + y and hi = 
xy, we can solve Eqs. (fl4|) and (fT5l) to obtain the exact 
number of spanning trees for the (x, y)-flowers: 

(x + y) n -(x + y-l)n-l 

N ST {n) = g n = {xy) t*+y-^ 

(x + B -2)( J: + B)"+(3! + B-l)n-( 3 i + ]< -2) 

x (x + y) {x+v-ir 2 (16) 
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FIG. 5: (Color online) The entropies of spanning trees in the 
(x, 2/)-flowers, depending on the two parameters x and y. 



For x = 1 and y = 2, Eq. (|T5|) recovers the result pre- 
viously obtained in [33]; for x = 1 and y — 3, Eq. (fT6|) 
reduces to Eq. (fTUI) ; while for a; = 2 and y = 2, Eq. (|T5)) 
agrees with the result given by Eq. ([Tl?]) . The consistency 
confirms the validity of the formula for spanning trees in 
the (x, j/)-flowers. 

Since the size for the (x, y)-flowers is N n — (x + y — 
2)/{x + y — l)(x + y) n + (x + y)/(x + y — 1), it is evident 
from Eq. (flU)) that NsT(n) scales exponentially with N n . 
This permits to get the entropies of spanning trees for 
the (x, y)-flowers as follows: 



E, 



(x : y)— flowers 



= lim 



hiNsT(ri) 



In xy + (x + y — 2) In (x + y) 
(x + y-l)(x + y-2) 



(17) 



It is worthwhile to mention that the above technique 
and process for finding spanning trees in the (x, j/)-flowers 
by evaluating determinants can be adapted to other self- 
similar networks. In order to show the universality, in the 
Appendix we apply the approach to determine the num- 
ber of spanning trees and their entropies for a small- world 
network (56j with an exponential deg ree distribution and 
a scale-free Apollonian network [57l l58j. and obtain ex- 
plicit solutions to the corresponding quantities. 



D. Analysis of the results 

The entropy for spanning trees is a very important 
quantity characterizing the network structure. Thus, it is 
interesting to show its dependence on the network struc- 
ture. Since, for the whole family of the (x, y)-flowers, its 
architecture is controlled by the two parameters x and 
y, it is expected that the entropy is dependent on x and 
y. In Fig. [SJ we show how the entropy varies with x and 
y. The intrinsic dependence relation of entropy on x and 
y is somewhat complex. Here we are mainly concerned 



about how the entropy -EY x ,y) -flowers is reliant on x and y 
when the sum x + y is fixed with x + y > 4. The main rea- 
son why we are interested in this case lies in that for this 
case it is meaningful to compare the entropy. As shown 
in section UH for the (x, y)-flowers with fixed x + y, they 
have an identical average degree. In addition, for fixed 
x + y, the networks have the same degree sequence and 
thus the same degree distribution. Particularly, for x > 1 
with fixed x + y, the networks are fractal having the frac- 
tal dimension monotonously decreasing in x. 

Figure [5] shows that when x + y is fixed, E/ Xiy \Q 0weIS 
is determined by the difference between x and y: the 
larger the difference, the smaller the entropy. Thus, in 
the case of fixed x + y, £'( x ,y)-flowera is an increasing func- 
tion of parameter x. For example, as shown above, when 
x + y = 4, the entropy for the (l,3)-flower is less than 
that for the (2, 2) -flower. Thus, we can conclude that the 
power-law degree distribution alone can not determine 
the entropies for spanning trees in scale-free networks. 
In fact, other than scale-free networks, the degree distri- 
butions of other networks do not suffice to characterize 
the number of spanning trees either. For instance, both 
the small- world network j56[ and the Koch network fam- 
ily [H, [6(j , to be addressed in the Appendix, have the 
same average degree and the same entropy of spanning 
trees, but their degree distributions are disparate: the 
former is exponential while the latter is scale-free. 

As pointed out above, degree distribution alone can- 
not determine the number of spanning trees. Then an 
important question arises: which structural property has 
the decisive effect on the number of spanning tress in 
the {x, y)-flower family with an identical value of x + yl 
Since, for the whole family of [x, y)-flowers with the same 
x + y, they have the same average node degree, the same 
degree sequence and thus the same degree distribution, 
and the same (zero) clustering coefficient, we argue that 
the difference in the number of spanning trees seen in 
Fig. [5] is attributed solely to the fractality. Specifically, 
the fractal property of networks significantly increases 
the number of spanning trees, which is related to the 
fractal dimension: the larger the fractal dimension, the 
less the number of spanning trees. This is obvious from 
Fig. [5] For those (x, y)-flowers with the same value of 
x + y, when x grows from 1 to its maximum, their fractal 
dimension df decreases from oo to its minimum, while the 
entropies of spanning trees increase from their minimum 
to maximum. Particularly, when x — 1, the network 
is non-fractal, it has the minimum number of spanning 
trees. 

Before closing this section, we discuss the case of differ- 
ent values of x + y. Figure [5] shows that for two network 
classes with different values of x+y, the class with smaller 
x + y has a larger entropy than that of the other class. 
This can be easily understood. Since for two networks 
with distinct values of x + y, their average node degrees 
are also disparate: the network with smaller x + y has 
a larger average node degree than that the other. Thus, 
the edge density of the former is greater than that of the 



latter, which is responsible for the difference of entropies 
for spanning trees in the two networks. 



IV. CONCLUSIONS 

In this paper, we have developed a rather generic tech- 
nique for evaluating determinants used to compute the 
numbers of spanning trees in self-similar networks. Our 
method is grounded on the recursive relations for the de- 
terminants of submatrices of the Laplacian matrix corre- 
sponding to the network at two consecutive generations, 
which arc established based on the structure of the net- 
work. An advantage of our approach is that it avoids 
the laborious computation of a determinant and reduces 
the computational complexity that a standard method 
needs. We have illustrated our technique by applying it 
to the family of (x, y)-flowers, which exhibits rich and 
special structural properties, as well as some remarkable 
features of many real systems. Moreover, to demonstrate 
the universality of our approach, we have applied it to 
some other networks with self-similarity. We obtained 
explicit formulas for the numbers of spanning trees in 
different networks. Our results show that degree distri- 
bution alone does not suffice to characterize the number 
of spanning trees in a network, and that fractal dimension 
has a predominant influence on the number of spanning 
trees in fractal scale-free networks with the same average 
degree. 
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Appendix: Spanning trees in other self-similar 
networks 



To illustrate the university of our proposed method, 
we apply it to calculate the number of spanning trees 
in some other self-similar networks with different degree 
distributions, exponential or scale-free. 



1. Small- world exponential network 

We first discuss a self-similar small- world network with 
an exponential form of degree distribution [56| , which is 
observed from some real-life systems 6l|. The small- 
world network is constructed iteratively (56[, denoted 
by S n after n (n > 0) iterations. The network starts 
from a triangle So, with three nodes called initial nodes. 
For n > 1, S n is obtained from S n —i: for each existing 
node in S n -i, two new nodes are generated, which and 




FIG. 6: (Color online) An illustration for the network S2. 



their mother node together form a new triangle. Figure[6] 
shows the structure of S2. It is easy to derive that the 
total number of nodes in S n is N n = 3™ +1 . The resultant 
network has a degree distribution decaying exponentially 
with the degree. In addition, it is small-world with the 
average distance growing logarithmically with its size. 

We now begin to evaluate the number of spanning trees 
in S n . In the case without any confusion, we use F n to 
denote the Laplacian matrix for S n , and use G n to denote 
a submatrix of F n , which is obtained by removing from 
F n the row and column corresponding to one of three 
initial nodes of S n . We have 
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(A.l) 



(A.2) 



where d n is the degree of an initial node in S n , and u n 
is a row vector representing the interaction between an 
initial node and the other non-initial nodes. 

According to the Matrix Tree Theorem, the number of 
spanning trees in S n is Nsrin) = det(G„). Using a sim- 
ilar iterative process as that in the main text, Eq. (|A.2[) 
leads to the following relation: det(G„) = 3(det(G„_i)) 3 . 
Considering the initial condition of dct(Go) = 3, we ob- 
tain the exact solution to the number of spanning trees 
in the small-world network S n as 



N ST (n) = 3< 3 " + - 1 )/ 2 
and the entropy of spanning trees as 



Eq = lim 



lnJVsxW = 1 
AT„ 2 n 



0.5493. 



(A.3) 



(A.4) 
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FIG. 7: The classic Apollonian packing of disks and its cor- 
responding Apollonian network. 

The above computational process for S„, can also be 
applied to the scale-free Koch networks [HH, [6(| with the 
same average node degree 3 as that of S n . Interestingly, 
the spanning tree entropy for the Koch networks is iden- 
tical to S n despite the fact that they have totally different 
types of degree distributions. 

2. Small-world scale-free Apollonian network 

We apply our technique to determine the spanning 
trees in the Apollonian network with both small-world 
and scale-free properties [13, H^|. The Apollonian net- 
work is derived from the classic Apollonian packing. This 
well-known packing problem [62[ starts with three mu- 
tually touching disks and their interstice is filled by a 



disk that touches all the three initial disks, forming three 
smaller interstices to be filled, as shown in the left panel 
of Fig. [7] Let each disk denote a node located at the 
center of the corresponding disk, linking the centers of 
touching disks by lines leads to the Apollonian network, 
as displayed in the right panel of Fig. [7] The Apollonian 
network can also be constructed iteratively and has the 
small- world effects [58|. After n iterations, the size of the 
Apollonian network is N n — 3 2 t " 5 . 

Making use of our technique for evaluating the cor- 
responding determinant, we can determine the number 
of spanning trees in the Apollonian network. Since the 
computational process is complicated and the relation 
equations are lengthy, we only provide the main results. 
The exact number of spanning trees in the Apollonian 
network after n iterations is 

N ST (n) = ~(3" + 5»)2 3 (3"-6n+3)/4 5 (3»-2„-l)/4 (A 5) 

Consequently, the entropy for spanning trees in the Apol- 
lonian network is 

£a P o1= lim ln y™) = ~ (In 3 + In 5) re 1.3540, 

(A.6) 

which is smaller than 1.9902, the entropy for spanning 
trees in three-dimensional lattices [7| having the same 
mean degree 6 as this Apollonian network. 
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